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Abstract The inclusive production cross sections of the strange vector mesons K*", K , 
and (f) have boon measured in interactions of 920 GoV protons with C, Ti, and W targets with 
the HERA-B detector at the HERA storage ring. Differential cross sections as a function of 
rapidity and transverse momentum have been measured in the central rapidity region and for 
transverse momenta up to pr = 3.5 GeV/c. The atomic number dependence is parametrised 
as CTpA = cTpN * , where apN is the proton- nucleon cross section. Within the phase space 
accessible, a{K*°) = 0.86 ± 0.03, a{K*°) = 0.87 ± 0.03, and a{(j>) = 0.96 ± 0.02. The total 
proton-nuclcon cross sections, determined by extrapolating the differential measurements 
to full phase space, arc a^jv^^.o = (5.06 ± 0.54) mb, a"pjv^i?*o = (4.02 ± 0.45) mb, and 
apN^tj) = (1.17 ± 0.11) mb. For all resonances the Cronin effect is observed; compared to 
the measurements of Cronin et al. for mesons, the measured values of a for 4> mesons 
coincide with those of K'^ mesons for all transverse momenta, while the enhancement for 

K*° /K*^ mesons is smaller. 

Key words. Proton nucleus collision, hadronic interaction, vector meson production, Cronin 

effect, atomic number dependence 



1 Introduction 

The search for a deconfined state of quarks and gluons in relativistic heavy ion colhsions is a 
subject which has received much attention in recent years [1,2]. It is a challenging task to find 
a signal that enables the discrimination of the short-lived locally deconfined state of quarks and 
gluons from a gas of confined hadrons. An enhancement of the strange quark fraction has been 
argued to appear in the quark gluon plasma [3] ; hence one expects an increase of strange hadron 
production in these processes. 

The rate of strange particle production depends on the tempcratTirc at which the hadroni- 
sation step takes place, and on the chemical potential. Hadronic interactions can influence the 
measurable strange particle rate in a two-fold way [4, 5] : they can induce transitions to other 
states, and, in the case of hadrons with a short lifetime, rescattering of the final state particles 
can take place, affecting the reconstructed invariant mass signal [6,7]. Since the strength of these 
final state interactions depends strongly on the lifetime, A'*'^ and (p mesons - where the lifetimes 
differ by a factor of about 10 - have been proposed as a versatile tool to study these effects in 
relativistic heavy ion reactions. 

Consequently, several new measurements of strange vector meson production have been pub- 
lished recently by RHIC experiments. Both for K*^ [8,9] and (p [10-16] mesons the influence of 
the very dense nuclear matter on the transverse momentum distributions and the line shape of 
the resonances have been studied in detail. In contrast to this, only limited information on (f) 
production in proton-nucleus interactions exists up to now. The ^d-dependence of (f) production 
in proton- nucleus interactions is studied in [17], while data in dAu collisions are available from 
the STAR and PHENIX collaborations [14,18,19]. No data are available on K*° production. 

Moreover, the study of the pr dependence of the inclusive production rate of K*'^ and (f) 
mesons in proton-nucleus (pA) interactions is of interest in its own right. Cronin et al. [20,21] 
observed a nuclear enhancement of high-pT particles in pA as compared to pp interactions for 
stable particles. In this paper we report the observation of the "Cronin effect" for the resonances 
K**^, K , and 0. The effect has been ascribed by many models to a broadening of the intrinsic 
transverse momentum of the primary partons in a hard scattering process [22] . The free param- 
eters are fixed by fitting the predictions to hard pp interactions. The high pT meson spectra for 
proton-nucleus interactions are reasonably well described by this ansatz, which however fails to 
reproduce the measured highpr baryon spectra [23]. On the other hand, the differences between 
meson- and baryon-production in proton-nucleus reactions can be explained by a final-state re- 
combination model [24]. Recently B.Z. Kopeliovich et al. [25] have formulated a parameter free 
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(a) (b) 

Figure 1. Invariant mass distributions for a) K^n^ and b) K^K~ mass combinations on the carbon target 
after application of the cuts described in the text. The distributions were fit by a Gaussian-smeared relativistic 
Breit-Wigner function and a background term (see text). 

description of the process. It explains the Cronin effect as a result of a convolution of the initial 
parton distribution affected by initial state interactions, followed by the hard scattering process 
and the fragmentation of the scattered partons. From this ansatz, it follows that the steeper the 
px dependence of the cross section, the larger the Cronin effect should be, i.e. its strength should 
depend on the flavor composition of the scattering partons via its structure- and fragmentation 
functions dependence. Moreover, a weak energy dependence of the Cronin effect is expected [25] 
in agreement with the observation of [21] for high-p^ mesons produced in proton-nucleus 
interactions in the energy range 200 GeV to 400 GeV. 

None of these models explicitly address the issue of the dependence of the Cronin effect 
on the flavor or mass of the produced hadron. The measurement presented here should moti- 
vate and aid further theoretical efforts to understand the Cronin effect and thereby lead to a 
deeper understanding of nuclear effects in pA collisions, particularly in the area of strangeness 
production. 



2 Data acquisition and event reconstruction 

The data were collected with the HERA-B detector operated at the proton ring of the HERA 
storage ring complex. The pN centcr-of-mass energy was i/i = 41.6 GeV. The K*^/K* data 
cover the rapidity interval —0.8 < y < 0.3 (in the centre-of-mass system) and a transverse 
momentum range of < < 12GeV^/c^. The (f) results were obtained for —0.7 < y < 0.25 
and 0.3 < Pt < 12GeV2/c2. Here the lack of data for pf, < O.SGeVVc^ is due to the small 
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opening angle of the (p decay in the lab system that makes the decay products escape through 

the b cam-pipe. 

The HERA-B detector is a forward magnetic spectrometer with a large acceptance for 
particles produced in the central region. The most important components of the detector for 
this analysis are a wire target [26] inserted into the halo of the 920 GeV proton beam, a vertex 
detector [27] followed by a magnetic dipole field of 2.13 Tm, a tracking system [28] using drift 
tube chambers, a Ring Imaging Cherenkov detector (RICH) [29] which enables the identification 
of mesons above a threshold of 9.6GeV/c, and an electromagnetic calorimeter [30]. A more 
detailed overview of the HERA-B detector can be found elsewhere [31]. 

The analysis is based on data recorded in the run period 2002/03. From a total of 200 • 10^ 
events collected with three different target materials a sub-sample of 132 • 10^ events taken in 
a 2-weeks period in December 2002 was selected to minimise systematic uncertainties. About 
58 ■ 10^ events were recorded with a carbon (C) wire target, 21 • 10^ events with a titanium 
(Ti) target, and 53 • 10^ events with a tungsten (W) target. As the data were recorded at a 
moderate interaction rate of 1.5 x 10^ s~^, compared to an average rate of filled proton bunches 
of 8.5 X 10^ s~^ at the target, only about 18% of all bunch crossings with the target produce 
an inelastic interaction. To reject empty bunch crossings an interaction trigger accepted only 
those events, where either the number of hits in the RICH was > 20 (for a (3 = 1 particle 32 
photons are detected on average) , or a cluster with an energy > 1 GeV was detected in the inner 
part of the electromagnetic calorimeter. The efficiency of the trigger was larger than 99% [32]. 
In parallel to the interaction trigger, events were selected in a random manner, regardless of 
whether or not they contained an inelastic interaction. These random events were used to study 
possible biases of the interaction trigger [32]. 

The events were reconstructed with the standard HERA-B analysis package [33], based on 
information from the vertex detector and the tracking system. Charged particles were identified 
by matching the tracks with rings reconstructed in the RICH. Using the track momentum 
information, the likelihood probability for charged kaons, protons, light particles (e, /x, tt^), and 
a background hypothesis were calculated [29]. 

The integrated luminosity was determined from the interaction rate at the target monitored 
by various devices and from the total inelastic cross section [32] . 

3 Data analysis 

3.1 K"^ meson identification and K*^ (^) optimisation 

To identify K*^/K*^ and (p, the decay channels K*^ K+ir" , Z*° ^ K-7T+ , and (p ^ K+K' 
were investigated. For all events with at least two oppositely charged tracks the invariant mass 
rriinv of the pair was separately calculated for the K~^7r^ , K^tt^, and K^K~ mass assumption. 
For the further analysis mass combinations were accepted if 

0.63GeV/c2 < rriinv < l-OGeV/c^ 

for the K^TT^, and 

0.988 GeV/c^ < rriinv < 1.198 GeV/c^ 

for the K~^K~ combinations, respectively. In order to reduce the combinatoric background, 
mesons were identified by the RICH. For the K*^{K* ) candidates, the kaon identification 
was demanded for the positive (negative) track, while for the (p candidates kaon identification 
was applied for both tracks. Only K"^ candidates with a momentum px > 9.7GeV/c were 
considered, slightly above the kaon threshold of the RICH. While for the determination of the 
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(j) signal a soft likelihood probability cut of Ck > 0.3 efficiently suppressed the background, a 
harder cut of Ck ^ 0.95 was required for the K*^ to reduce the combinatoric background. The 
non-existence of mass reflections has been checked by experimental data and by the Monte-Carlo 
(MC) simulation described in chapter [331 

For momenta pK < 40 GeV/c the efficiency for identification of mesons is high, amounting 
to about 90 % [29]. For higher momenta the discrimination power drops and the misidentification 
probability increases, since for such momenta the Cherenkov angles for tt and K mesons become 
indistinguishable. The K meson identification efficiency as a function of K momentum has been 
determined by MC. It was checked with (p K^K~ data [34] by identifying either one or both 
kaons and reconstructing the (j) resonance. The difference in efficiency between data and MC of 
5% per identified kaon was corrected for, and a residual systematic error of 2% was estimated 
(see chapter [3^ . A large fraction of the produced K*^ and </> candidates were rejected by the 
momentum cut pK > 9.7GeV/c. The total efficiency for K mesons from K*^ and (j) resonance 
decays to survive the RICH identification and the momentum cut was about 15 %. 

The significance of the signal was improved further by applying soft track quality cuts. In 
case of the vertex detector at least 6 hits per track were demanded. In the tracking system 
a cut of at least 15 hits per track was used to enhance the significance of the K*^ ((p) signal 
without reducing the efficiency. On average, a track comprised 12 hits in the vertex detector and 
40 hits in the tracking system, respectively. The K*^ signal significance was improved further 
by applying a cut of Pt^ > 2GeV/c to pion track candidates [35]. The cuts were optimised for 
significance using signal MC and real data background events. 



3.2 Signal determination 

Typical invariant mass distributions for the K^Tr~ and K^K^ combinations are shown in Fig. 
[U after application of the cuts mentioned above. To facilitate the extraction of the signal, the 
resonance was described by a relativistic P-wave Breit-Wigner [36]: 

BW{m) = (1) 

(mg — m^Y -h (mo • l \m)Y 

rM = r„. (A) (2) 

with 

Fq = r{mQ) the natural width 
g(m)|^.o 



(m^ — {niK + mT^Y){rn? — {mx — rn.,^)'^) 



2m 



2 

go = qimo) 
I = 1. 

q{ra) is the momentum of the decay products in the rest frame of the mother particle, and / is 
the angular momentum transferred. To take into account effects of the detector resolution the 
Breit-Wigner distribution was numerically folded with a Gaussian of width am in the fits of 
the measured mass spectra. It has been checked that the mass resolutions in experimental data 
and MC are consistent within errors. The mass mo of the resonance and the width am were 
free parameters of the fit, while the natural width Fq of the resonance was fixed to its PDG 
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Table 1. Number of inelastic events (Nia), number of if*", K*'^ , and tj) candidates (Nvm) after analysis cuts, 
reconstructed mass (mrec), and experimental mass resolution (um) for each target material. Errors are statistical 
only. 







Nvm 


rur^c [MeV/c^] 


a™ [McV/c^] 


pA K*°X 


C 


58- 


10® 


187574 ± 1726 


890.5 ±0.2 


5.1 ± 0.2 


Ti 


21 • 


10*^ 


83231 ± 1237 


890.7 ±0.4 


5.0 ±0.3 


W 


53 • 


lO** 


262356 ± 2420 


890.8 ±0.2 


5.3 ±0.3 








pA - 






C 


58- 


10® 


149721 ± 1622 


890.8 ±0.2 


4.6 ±0.1 


Ti 


21 ■ 


10® 


65823 ± 1147 


890.7 ±0.4 


5.1 ± 0.4 


W 


53- 


10® 


195373 ± 2225 


891.5 ±0.3 


4.8 ±0.2 








pA 






c 


58- 


10® 


15379 ± 204 


1019.2 ±0.1 


2.6 ±0.1 


Ti 


21 • 


10® 


7580 ± 142 


1019.2 ±0.1 


2.7 ± 0.1 


W 


53- 


10® 


26971 ± 295 


1019.3 ±0.1 


2.6 ±0.1 



Table 2. Summary of systematic errors for the integrated and total production cross section measurements. See 
text for a more detailed explanation. Numbers are grouped target material wise. The last two errors are due to 
the phase space extrapolation and have been applied only to the total cross sections. Numbers are given in %. 







C 






Ti 






W 










<^ 














MC statistics 


1.0 


1.3 


1.6 


1.0 


1.3 


1.1 


1.0 


1.3 


0.7 


Bg. parametrisation 


2.2 


1.9 




1.7 


2.3 




2.5 


2.3 




Cut variation 


2.8 


2.8 


6.2 


2.8 


2.8 


7.0 


2.8 


2.8 


7.1 


BR 


0.1 


0.1 


1.22 


0.1 


0.1 


1.22 


0.1 


0.1 


1.22 


Luminosity 




4.8 






5.0 






3.8 




Luminosity scaling 




2.0 






2.0 






2.0 




Track efficiency 




3.0 






3.0 






3.0 




p|n-extr. 






6.6 






8.8 






8.2 


y-exti. 


2.6 


4.1 


6.3 


2.6 


4.5 


6.0 


2.6 


4.5 


4.9 



value [37]. The number of resonance decays was estimated by integrating the fitted signal in 
the interval ±3/0 around the resonance mass ttiq. A variation of the fit interval by Fq changed 
the yield up to 5%. The same was observed for the MC used for acceptance correction, thus 
compensating the effect seen in data. 

The background was parametrised by expressions which have been applied successfully in 
other experiments. For the K*^/K* the parametrisation of [38] was adopted: 

BGK,o{m) = po ■ (^Y" ■ e-P'" " ^'3^', (3) 
Vm/ 

where pQ . . .p^ are free parameters and q is the momentum transfer introduced above. For the 
(j) meson the background was described by a parametrisation proposed in [39]: 

BG^m) = po ■ mZ ■ e'^^"^ " p^^' , (4) 

where rrirei = m — 2mK and po • • - Ps again are free parameters of the fit. As shown in Fig. [U 
these combinations of signal and background functions describe the data properly. 

It has been shown [34,35] that both the fitted masses and the detector resolution am derived 
from the data were stable for the whole data taking period and independent of the target 
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material. Also the resonance yields, normalized to the integrated luminosity, were stable during 
the whole data taking period. Table [T] summarises the i^*^, K* and (p statistics obtained after 
cuts, the reconstructed mass positions, and the experimental mass resolutions. 

The bin size for the p^- and rapidity (y-) distributions were chosen to be broader than the 
resolution determined from MC simulation. From MC we find that in each bin more than 80% 
of the reconstructed vector mesons were generated in the same bin. 



3.3 Acceptance correction 

Acceptance correction functions were determined from MC events using FRITIOF 7.02 [40]. The 
detector simulation was based on GEANT 3.21 [41] and includes the measured hit resolutions, 
mapping of inefficient channels, and electronic noise. The simulated events were processed by the 
same reconstruction chain as the data. Since the FRITIOF generator uses an approximation for 
the generation of hard parton-parton scattering, the MC p'^ distributions show a steeper slope 
than observed in the data. To make sure that the reconstruction efficiencies and the influence 
of the analysis cuts are correctly described by the simulation, the MC events were reweighted 
in p'^ according to the data for both the K*^ / K and the cj) meson. It was checked that after 
reweighting both the rapidity and azimuthal MC distributions match those of the data. 

The total efficiency for each p^- and y-bin was determined with the help of the reweighted 
MC events by comparing the number of reconstructed and produced K*^/(j) mesons in each bin. 
The number of reconstructed mesons was obtained similar to data by a fit to the reconstructed 
invariant mass spectra. Typical efficiencies vary between 5% and 30% [34,35]. Studies of sys- 
tematic uncertainties due to signal extraction and background estimation show that they are 
negligible compared to the statistical error in the corresponding interval. Note that for each 
rapidity bin the efficiencies are determined by integrating over the measured p'^-range and vice 
versa. 



3.4 Systematic errors 

Possible sources of systematic uncertainties have been studied in detail: Since especially the 
K*^/K* invariant mass distributions are affected by a large combinatorial background peaking 
close to the resonance signal, an alternative background parametrisation to ([3]) was used, which 
was proven to be successful in studies of K*^ production in e"'"e~-collisions [42]: 



BG{m) = po 



1 



1 



m+ pi 



P2 



[l+psm), (5) 



where po • • - Ps are free parameters. The signal yields obtained in most pj,- and y-intervals agree 
for the two analyses within their statistical uncertainty. The absolute difference of the yields, de- 
rived using the two background parametrisations, is considered as a systematic error. It amounts 
to about 2% for the integrated cross sections. For the cp meson a corresponding study [34] pro- 
vided much smaller deviations, and their contribution to the systematic error can be safely 
neglected. 

The variation of the analysis cuts, including the uncertainties in the kaon identification 
efficiency (2.0% per identified particle), resulted in a systematic uncertainty of 2.8% for the 
K*^/K*^ and 6.2-7.1% for the cj) analysis (depending on the target material). The systematic 
error of the luminosity consists of an overall scaling error of 2% and a target material dependent 
part of 3.8-5.0%, and the systematic error of the branching ratios [37] amounts to 0.1% and 
1.22% for K*^ K+ir-fK K"7r+ and (p K+K', respectively. A systematic error on 
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Figure 2. Measured inclusive differential cross sections da/dy as a function of rapidity for a) if*" and b) K*" 
meson production in the range < pj- < 12 GeV^/c'^. Horizontal bars indicate the widths of the individual bins, 
while vertical errors bars represent the uncertainties due to statistics. Note that for most bins the error is smaller 
than the size of the symbols. The dashed lines show the (scaled) results of the FRITIOF MC. 



the track efficiency of 1.5% per track was estimated using 
of all systematic errors is given in Table [21 



decays [43]. A summary 



4 Experimental results 

4.1 Differential cross sections 



Due to limited statistics, we did not calculate the double-differential cross section, but rather 
determined the one-dimensional rates dN/dp^ and dN/dy, which were integrated over the mea- 
sured y and p'^ range, respectively. Given the corresponding one-dimensional efficiencies 
eacc(Pj') and eacciv), the differential cross sections can be written as: 



da _ 1 dN 1 
dx BR ■ C dx e„rr( 



(6) 



where x = BR and C denote the branching ratio of the decay and the integrated lumi- 

nosity of the data set, respectively. The acceptance boundaries of the measurement in rapidity 
and transverse momentum squared are given in Table [H 

The measured cross sections for the K*^ and the K*^ mesons as a function of rapidity are 
shown in Fig. [2l The cross sections are rather flat as expected in the central rapidity region [44] . 
The K*^ cross sections are suppressed in comparison to the K*^ cross sections by about 20%, 
in agreement with the results of the FRITIOF MC. The rapidity distributions for (j) mesons 
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produced on C, Ti, and W targets are shown in Fig. [H Again the cross section in the central 
region is flat. 

Transverse momentum spectra are shown in Fig. [3] for the K*^ and K* . The spectra deviate 
from a simple exponential shape at large transverse momenta, which can be traced back to hard 
parton-parton scattering [45]. Common exponential parametrisations like 

dpi ' ' ^' 

or 



da 



Cse-^™^ mT = Jml+pl (8) 



dp 



,2 



only fit the data for < 1 GeV^/c^. 

On the contrary, as is shown in Fig. [3l a power-law parametrisation, 

describes the spectra over the full p^ range of the measurements. The p|i-spectrum of K*^ 
mesons, which have no valence quarks in common with nucleons, is steeper than the K*^ meson 
spectrum. Similar observations have been made in hard pp collions for K~ mesons in comparison 
to mesons [45,46]. In the latter case this observation has been traced back to the fact that K~ 
mesons are, for the kinematic condition of the experiment, dominantly gluon fragments [45,47]. 
QCD calculations [48] show that for a fixed value of xt = 2 ptI \fs the relative contribution of 
gluon fragments to K~ {K ) production decreases while the sea quark share is enhanced when 
the CMS energy is increased. 

The same parametrisation ([9]) successfully fits the transverse momentum distributions for 
the (j), shown in Fig. [5l The suppression of (f) mesons as compared to K mesons by a factor of 
3 to 4 is compatible with the known strangeness suppression in fragmentation processes [42]. 
No simple qualitative explanation exists for the observation that the K px spectrum is steeper 
than that of (p meson. 

The fit parameters for all particles are collected in Table [3l Note that the of the tungsten 
fits is considerably worse than those of the other target materials. This is likely caused by an 
enhancement of background due to the large multiplicity of the pW interactions. The improves 
however considerably when leaving out the bin of largest p^ in the fits. The differential cross 
sections da/dy and da/dp^ for K*^ , K*^ , and (j) production are summarised in Tables [6] and [71 



4.2 Integrated inclusive and total cross sections 

The integrated cross sections in the acceptance region of the HERA-B detector are calculated 
by fitting the shape of the MC distributions to the differential cross sections da/dy. The results 
of the fits are shown as dashed lines in Figs. [2]and|H The integrated cross sections ams in the 
phase space accessible to the HERA-B experiment are listed in Table [H Note that integrating 
the p1 differential distributions yields consistent cross sections. 

The good agreement of the K*^ and K*^ transverse momentum spectra with the power-law 
parametrisation ([9]) for small pT suggests that this is true also for the (p differential cross sections, 
where only data with p^ > 0.3GeV^/c^ are available. Hence we have used this parametrisation 
to extrapolate the cp differential cross section to < < 0.3 GeV^/c^ when calculating the total 
cross section. The extrapolation factors are summarised in Table [H The systematic uncertainty 
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Figure 3. Measured inclusive differential cross sections da/dp^ as a function of for a) JsT*" and b) K*" 
production. The vertical error bars reflect statistical errors only. Fits to the parametrisation (|9} are superimposed. 
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Figure 4. Measured inclusive differential cross section da/dy as a function of rapidity for (p meson production in 
the range 0.3GeV^/c^ < Pt < 12GeV^/c^. The vertical error bars reflect the uncertainties due to statistics. The 
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Figure 5. Measured inclusive differential cross sections da/dp^ as a function of for (j) meson production. The 
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Table 3. Results of the power-law fits ^ for the differential cross sections da/dp^. Fits were performed within 
the full pl^-range of the measurements and in the rapidity interval accessible by the experiment (—0.8 < y < 0.3 



for K'°/K , -0.7 <y < 0.25 for 
the fit parameters. 



The last column gives the mean pr of the cross sections, calculated from 





Ci [mb/GeV^] 




pI [GeVVc'] 


xVn.d.f. 


(pt) [GeV/c] 


pA K*°X 


C 


28.5 ± 1.1 


4.17 ± 0.12 


1.36 ±0.08 


9.1/8 


0.66 ±0.07 


Ti 


87.1 ±5.2 


3.89 ±0.14 


1.30 ±0.11 


8.3/8 


0.69 ±0.09 


W 


268.5 ± 10.9 


3.75 ±0.10 


1.32 ±0.08 


24.5/8 


0.72 ±0.07 


pA K*°X 


C 


24.4 ± 1.0 


4.67 ±0.16 


1.52 ±0.11 


14.2/8 


0.64 ±0.07 


Ti 


71.9 ±4.4 


4.39 ±0.21 


1.51 ±0.14 


13.0/8 


0.67 ±0.11 


W 


230.5 ± 10.2 


4.51 ± 0.16 


1.68 ±0.12 


36.9/8 


0.69 ±0.08 


pA^ (j>X 


C 


7.6 ±0.8 


4.04 ± 0.15 


1.26 ±0.13 


12.6/12 


0.66 ±0.09 


Ti 


31.9 ±5.0 


3.81 ± 0.21 


1.15 ±0.17 


10.5/12 


0.66 ±0.14 


W 


73.9 ±6.1 


4.20 ±0.14 


1.65 ±0.14 


28.4/12 


0.72 ± 0.09 



of the extrapolations is estimated by using parametrisation d?]) as an alternative and taking 
the difference between the two approaches. Leaving out the bin of largest in the power-law 
fit increases the extrapolated (p cross section by about 6%, which is fully compatible with this 
systematic uncertainty. 

Taking into account the limited rapidity coverage of the detector, the integrated cross sections 
are then extrapolated to full phase space using target-specific correction factors Aams- These 
factors are derived from the MC simulators FRITIOF [40] and HIJING [49], and the average 
of both was taken for extrapolation. The systematic error was estimated by taking half the 
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Table 4. Extrapolated total cross sections (JpA, and integrated inclusive differential production cross sections avis 
within the accessible phase space of the measurement. The errors denote the statistical error and the combined 
systematic error, respectively. The rapidity and coverage of the measurements are given in columns 3 and 
4. For the (j> meson, R'^f'^ denotes the factor used for extrapolation of the integrated cross sections to = 0. 
Aovis is the fraction of the total cross section within the rapidity coverage of the measurement, based on the 
FRITIOF [40] and HIJING [49] Monte-Carlo generators. Note that for these numbers the py-extrapolation for 
the 4> cross sections was applied beforehand. 





apA [mb] 


cr„is [mb] 


rapidity interval 


p% range 


[GeVVc''] 




Aay,s [%] 


pA K*°X 


C 


43.9 ±0.6 ±3.3 


12.1 ±0.2 ±0.9 


-0.8 - +0.3 


0.0 - 


- 12.0 


1 


27.4 ±0.7 


Ti 


141.2 ± 2.6 ± 10.6 


38.5 ±0.7 ±2.7 








1 


27.3 ±0.7 


W 


465.9 ±6.4 ±32.7 


127.5 ± 1.7 ±8.3 








1 


27.4 ±0.7 


pA K*°X 


C 


36.0 ±0.6 ±2.9 


10.0 ±0.2 ±0.7 


-0.8 - +0.3 


0.0 - 


- 12.0 


1 


27.7 ± 1.1 


Ti 


111.5 ±2.5 ±9.7 


31.1 ± 0.7± 2.3 








1 


27.9 ± 1.3 


W 


388.8 ±6.9 ±30.8 


107.6 ±1.9 ±7.0 








1 


27.7 ± 1.3 


pA^ (j)X 


C 


11.0 ± 0.2 ± 1.4 


1.60 ±0.03 ±0.14 


-0.7 - +0.25 


0.3 - 


- 12.0 


1.94 + 0.13 


28.3 ± 1.8 


Ti 


44.7 ± 1.0 ±6.0 


6.4 ±0.1 ±0.6 








1.97 ±0.17 


28.3 ± 1.7 


W 


135.1 ± 1.6 ± 17.0 


21.9 ±0.3 ±2.0 








1.72 ±0.14 


27.9 ± 1.4 



Table 5. Atomic mass number dependence of the integrated and total production cross sections, and extrapolated 
production cross sections in proton-nucleon reactions. The errors of a reflect the statistical errors and all target- 
material dependent systematic errors. For o-pjv, full errors due to statistics and systematics are given. Errors were 
added in quadrature. 





pTV ^ K'^^X 


pN ^ 


pX (j)X 




0.86 ±0.03 


0.87 ±0.03 


0.96 ±0.02 


a-pA 


0.87 ±0.03 


0.87 ±0.03 


0.91 ±0.02 


(jpN [mb] 


5.06 ±0.54 


4.02 ± 0.45 


1.17±0.11 



difference between the FRITIOF and the HIJING results. As quantified in Tabled about 28% 
of the total cross section is accessible to the experiment. Correcting the limited phase space 
coverage leads to the total inclusive cross sections 

T> extr 

■ w vis 

= -^^^ (10) 

VIS 

for K*^ ^ K*^ , and (p production in pC-, pTi-, and pW-interactions (collected in Tabled]). The 
systematic error of the cross section extrapolation in total amounts to about 2.6% for the K*^, 
4.5% for the K , and 10.0% for the (j) total production cross sections (Table [2|). 



4.3 Atomic mass number dependence of the integrated and total cross sections 

Fig. [7] shows the atomic mass number dependence of the integrated cross sections avis for the 
production of K*^, K* , and (f) mesons in pA-interactions. The atomic number dependence is 
well-described by a power-law 

CT^i5 = • (11) 

where avis,o is the visible proton-nucleon cross section. The results for the exponent Ovis are 
collected in Table [5l Comparing the values with those for inelastic pA cross sections, amei = 
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Figure 6. a) Comparison of K*'^ /K* and K*^ meson total production cross sections at various centre-of-mass 
energies [38,51-57]. The results of this analysis have been slightly displaced from each other for better visibility, b) 
same as a) but for (p meson production [17,38,54,58-63]. Results of the pBe measurements have been extrapolated 
to A = 1 using a = 0.91 ± 0.02 from [17]. 



0.71 ±0.01 [50], it follows that strange vector meson production shows a stronger 74-dependence 
than the total inelastic cross section. While the latter scales roughly with the cross-sectional 
area of the nucleus, the K*^ and (f) production cross sections show a tendency towards scaling 
with the volume. 

The same parametrisation can be used to describe the ^-dependence of the cross sections 
extrapolated to the full phase space, 

apA = CTpN-A''r>A^ (12) 

is an estimate of the total cross section for strange vector meson production in proton- 
nucleon reactions. The fitted values for apA and apN are included in Table El and the results for 
(TpAT are compared to those measured at other energies in pp and pA reactions (Fig. [6]) . Recently, 
the NA60 collaboration measured the production of (p mesons in collisions of 400 GeV protons 
with various nuclei [17]. Fitting their data to eq. [T2l they obtain a = 0.91 it 0.02, which is in 
good agreement with our measurements. Consequently, the pBe data displayed in Fig. [6]b) were 
rescaled using this result. 



4.4 Cronin effect 

The A-dependence of the differential cross sections da/dy and da/dp^ themselves is also of 
interest. Fig. [8] demonstrates that, for the differential cross sections integrated over p^, the 
exponent a is in good approximation independent of rapidity in the range of the measurement. 
Note that mainly soft interactions with < IGeV^/c^ contribute to the cross section. 

In contrast, a depends on the transverse momentum. With increasing p^, a constantly in- 
creases, as shown in Fig. O Cronin et al. [21] have observed this effect for the production of 
stable particles (vr, K, p). In this analysis we report an observation of the Cronin effect for the 
production of vector mesons containing strange quarks. Note however, that the exponent for 
K*^ and K* production is smaller than that for and K~ mesons. Only at the highest pT 
values do the distributions get close to each other, while the exponents a for (j) and mesons 
coincide for all transverse momenta. 

Qualitatively the observations shown in Fig. [9] are in agreement with the model of [25] which 
predicts a decrease of a with increasing CMS energy; this is indeed observed for the K*^ (K* ) 
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Figure 8. Measured values of ct{y) for a) K*^ jK and b) production. 



data collected at a beam energy of 920 GeV when compared to the (-f^ ) results of Cronin et 
al. measured at 400 GeV. On the other hand produced at a beam energy of 920 GeV coincides 
with that of mesons collected at 400 GeV; this might indicate a different flavor composition 



of the initially scattered partons fragmenting to (/> mesons as compared to i^* {K ) mesons. 



5 Conclusion 



The cross sections of K*^ ^ K*^ , and (j) meson production in the central rapidity region have 
been measured for the first time in pC-, pTi-, and pW-inter actions. The cross sections show 
a power-law dependence a oc A°' with an exponent of a « 0.8 — 1.1 depending both on the 
produced vector meson and the transverse momentum. The exponent a is found to be larger 
than the corresponding value for the total inelastic cross section. The measurements of a{p^) 
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Figure 9. Measured values of q(pt) for a) A"*", b) K and c) (j) production. The results of Cronin et al. [21] on 
the mesons are superimposed. Numbers of this analysis are summarised in Table [T] 



for the (p vector meson agree with the values determined by Cronin et al. [21] for mesons. 
For the K*^/K mesons, a(pj') is systematically smaller. 
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Table 6. Compilation of the differential production cross sections da/dy in the accessible phase space of the 
measurement. The errors reflect the statistical and full systematic errors, respectively. 





Ay 




da/dy [mb] 










c 


Ti 


W 


K 


-0.8 


-0.6 


12.5 ± 0.4 ± 0.9 


38.6 ± 2.2 ± 2.9 


128.3 ± 4.6 ± 8.6 




-0.6 


-0.5 


12.3 ± 0.5 ± 0.9 


A -i 1 n 1 on 

41.6 ± 2.2 ± 3.2 


131.6 ± 5.4 ± 9.0 




-0.5 


--0.4 


10.8 ± 0.4 ± 0.8 


34.0 ± 1.8 ± 2.5 


1 orr 1 ^o 

112.9 ± 3.7 ± 7.3 




n A 


— u.o 


Iz.z ± U.o ± u.y 


Oil. 1 ± L.I ± z.y 


izU.D ± 0.4 ±1.1 




-yj.o 


n 9 


11.1 ±0.3 ±0.8 


34.4 ± 1.7 ± 2.4 


114.3 ±3.0 ±7.1 




n 9 


— -U. 1 


10.6 ±0.3 ±0.7 


33.5 ± 1.4 ±2.3 


112.8 ±3.0 ±7.0 




n 1 

-U. i 


— u.u 


10.7 ±0.3 ±0.7 


34.6 ± 1.4 ±2.4 


115.2 ±2.9 ±7.1 




n n 
u.u 


-|-U.i 


11.1 ±0.3 ±0.8 


33.5 ± 1.5 ± 2.4 


121.2 ±3.3 ±7.6 




-hU.l 


J-O 9 


10.1 ±0.3 ±0.7 


32.4 ±1.5 ±2.4 


110.3 ± 3.6 ± 7.2 




-1-0 9 


i_n 


10.0 ±0.4 ±0.8 


37.0 ±2.3 ±2.8 


109.0 ± 3.9 ± 7.5 


K 


-0.8 


--0.6 


8.5 ± 0.4 ± 0.7 


28.7 ± 1.7 ± 2.3 


96.8 ± 4.1 ± 7.1 




-0.6 


--0.5 


9.9 ± 0.5 ± 0.8 


30.0 ± 2.1 ± 2.5 


114.0 ± 4.9 ± 8.6 




-0.5 


--0.4 


9.2 ± 0.4 ± 0.7 


31.5 ± 1.8 ± 2.5 


97.7 ± 4.3 ± 6.9 




-0.4 


- -0.3 


Q o 1 n /I 1 n ^ 
o.o ± U.4 zh U. / 


ZO.D ± i.D ± Z.i 


1 no 4-/1 4-7/1 
iUz.y zt 4.U 4= (.4 




-0.3 


- -0.2 


9.4 ± 0.3 ± 0.7 


31.1 ± 1.5 ± 2.3 


102.8 ± 3.4 ± 6.7 




-0.2 


- -0.1 


9.3 ± 0.3 ± 0.7 


31.5 ± 1.8 ± 2.4 


105.6 ± 3.5 ± 7.0 




-0.1 


- 0.0 


9.5 ± 0.3 ± 0.7 


28.4 ± 1.5 ± 2.1 


91.7 ± 3.3 ± 6.1 




0.0 


+0.1 


8.7 ± 0.3 ± 0.6 


26.0 ± 1.3 ± 1.9 


88.8 ± 3.4 ± 5.6 




-1-0.1 


- -1-0.2 


8.9 ± 0.3 ± 0.7 


25.1 ± 1.6 ± 2.0 


97.3 ± 3.2 ± 7.2 




-1-0.2 


- -1-0.3 


8.1 ± 0.4 ± 0.7 


25.9 ± 1.8 ± 2.2 


92.3 ± 4.0 ± 7.1 




-0.70 


- -0.55 


1.77 ± 0.10 ± 0.18 


7.2 ± 0.6 ± 0.8 


23.1 ± 0.9 ± 2.1 




-0.55 


- -0.40 
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i.ol zt U.UO zt U.i / 


70-1-0/1-1-07 
/ .U zt U.4 zt U. / 


9/1 0-1-0(^-1-9 1 
Z4.U zt U.D zt Z.l 




-0.40 


- -0.30 


1.80 ±0.06 ±0.21 


7.2 ±0.3 ±0.7 


24.7 ±0.7 ±2.2 




-0.30 


- -0.20 


1.60 ±0.05 ±0.14 


6.8 ±0.3 ±0.7 


23.0 ±0.5 ±2.0 




-0.20 


- -0.10 


1.71 ± 0.05 ± 0.14 


6.7 ±0.3 ±0.6 


21.7 ±0.5 ±1.9 




-0.10 


- 0.00 


1.64 ±0.05 ±0.14 


6.4 ±0.3 ±0.6 


22.7 ±0.6 ±2.0 




0.00 - 


- +0.W 


1.84 ±0.07 ±0.16 


7.3 ±0.4 ±0.7 


23.6 ±0.8 ±2.1 




-hO.lO 


- +0.25 


1.45 ±0.10 ±0.14 


6.1 ±0.5 ±0.6 


26.5 ±1.1 ±2.5 
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Table 7. Compilation of the differential production cross sections da/dp^ in the accessible phase space of the 
measurement. The errors reflect the statistical and full systematic errors, respectively. The right column lists the 
results of the A°' flts shown in Figs.|9l The errors given are the combined statistical and target-material dependent 
systematic errors, added in quadrature. 



A{p^) [GeVVc^] da/dp^T [mb/(GeV/c)^] 









C 


Ti 


w 


a 




0.0 


- 0.2 


22.2 ± 0.5 ± 1.5 


69.4 ± 4.3 ± 4.9 


224.3 ± 5.8 ± 14.0 


0.85 ±0.03 




0.2 


- 0.4 


12.2 ±0.4 ±0.8 


38.9 ±1.6 ±2.7 


119.4 ±4.2 ±7.4 


0.84 ±0.03 




0.4 


- 0.6 


7.4 ± 0.2 ± 0.5 


25.5 ± 1.3 ± 1.8 


76.0 ± 2.7 ± 4.9 


0.85 ± 0.03 




0.6 


- 0.8 


4.7 ± 0.2 ± 0.3 


14.5 ± 0.8 ± 1.1 


50.5 ± 2.0 ± 3.3 


0.87 ± 0.03 




0.8 


- 1.0 


3.5 ± 0.1 ± 0.3 


10.6 ± 0.7 ± 0.8 


37.9 ± 1.3 ± 2.5 


0.88 ± 0.03 




1.0 


- 1.5 


1.92 ± 0.04 ± 0.13 


6.4 ± 0.2 ± 0.4 


21.8 ± 0.5 ± 1.3 


0.89 ± 0.02 




1.5 


- 2.0 


0.91 ± 0.03 ± 0.07 


3.1 ± 0.2 ± 0.2 


11.9 ± 0.3 ± 0.8 


0.95 ± 0.03 




2.0 


- 3.0 


0.37 ± 0.01 ± 0.03 


1.45 ± 0.06 ± 0.11 


5.0 ± 0.1 ± 0.3 


0.95 ± 0.03 




3.0 


- 4.0 


0.157 ± 0.006 ± 0.013 


0.56 ± 0.04 ± 0.05 


2.20 ± 0.08 ± 0.16 


0.97 ± 0.04 




4.0 


- 7.0 


0.035 ± 0.001 ± 0.003 


0.142 ± 0.008 ± 0.011 


0.64 ± 0.02 ± 0.04 


1.08 ± 0.03 




7.0 


- 12.0 


0.0044 ± 0.0003 ± 0.0004 


0.021 ± 0.002 ± 0.002 


0.084 ± 0.004 ± 0.007 


1.08 ± 0.05 


K*° 


0.0 


- 0.2 


19.3 ± 0.5 ± 1.4 


59.2 ±3.4 ±4.3 


201.8 ±5.5 ±13.0 


0.86 ±0.03 




0.2 


- 0.4 


10.1 ±0.3 ±0.7 


31.0 ±1.5 ±2.2 


101.3 ±3.7 ±6.5 


0.85 ±0.03 




0.4 


- 0.6 


6.2 ±0.2 ±0.4 


19.1 ± 1.1 ± 1.4 


65.0 ±2.9 ±4.3 


0.86 ±0.03 




0.6 


- 0.8 


4.0 ±0.2 ±0.3 


13.9 ± 0.9 ± 1.1 


46.2 ±2.3 ±3.2 


0.90 ±0.04 




0.8 


- 1.0 


2.8 ±0.1 ±0.2 


10.1 ±0.7 ±0.8 


34.5 ± 1.7 ± 2.5 


0.93 ±0.04 




1.0 


- 1.5 


1.55 ±0.04 ±0.12 


5.1 ±0.3 ±0.4 


19.0 ±0.7 ±1.3 


0.92 ±0.03 




1.5 


- 2.0 


0.68 ±0.02 ±0.05 


2.6 ±0.1 ±0.2 


8.7 ±0.3 ±0.6 


0.93 ±0.03 




2.0 


- 3.0 


0.264 ±0.010 ±0.022 


0.88 ±0.05 ±0.07 


3.9 ±0.1 ±0.3 


1.00 ± 0.04 




3.0 


- 4.0 


0.102 ±0.006 ±0.011 


0.45 ± 0.03 ± 0.04 


1.78 ±0.08 ±0.16 


1.05 ±0.05 




4.0 


- 7.0 


0.024 ± 0.001 ± 0.002 


0.097 ± 0.008 ± 0.009 


0.42 ±0.02 ±0.03 


1.06 ±0.05 




7.0 


- 12.0 


0.0018 ± 0.0002 ± 0.0002 


0.010 ± 0.002 ± 0.001 


0.035 ± 0.003 ± 0.003 


1.07 ±0.07 




0.3 


- 0.4 


2.8 ±0.3 ±0.2 


9.4 ± 1.7 ± 1.2 


34.9 ±4.1 ±3.8 


0.92 ±0.07 




0.4 


- 0.5 


2.3 ±0.2 ±0.3 


9.2 ±1.3 ±1.1 


31.2 ±2.5 ±3.1 


0.95 ±0.06 




0.5 


- 0.7 


1.70 ±0.09 ±0.16 


7.0 ±0.6 ±0.7 


23.0 ± 1.1 ± 2.1 


0.95 ±0.03 




0.7 


- 0.9 


1.09 ±0.05 ±0.09 


4.5 ±0.3 ±0.5 


14.2 ±0.6 ±1.3 


0.94 ±0.03 




0.9 


- 1.1 


0.73 ±0.03 ±0.07 


2.8 ±0.2 ±0.3 


10.1 ±0.4 ±0.9 


0.96 ±0.03 




1.1 


- 1.3 


0.46 ±0.02 ±0.04 


2.0 ±0.1 ±0.2 


6.6 ±0.2 ±0.6 


0.97 ±0.03 




1.3 


- 1.5 


0.35 ±0.02 ±0.03 


1.57 ±0.09 ±0.16 


5.2 ±0.2 ±0.5 


0.98 ±0.03 




1.5 


- 1.8 


0.266 ±0.010 ±0.024 


1.02 ±0.06 ±0.10 


4.0 ±0.1 ±0.4 


1.00 ± 0.03 




1.8 


- 2.2 


0.163 ±0.006 ±0.015 


0.62 ± 0.03 ± 0.06 


2.59 ±0.07 ±0.23 


1.02 ±0.03 




2.2 


- 2.7 


0.098 ± 0.001 ± 0.009 


0.43 ± 0.02 ± 0.04 


1.57 ±0.04 ±0.14 


1.01 ±0.02 




2.7 


- 3.3 


0.058 ± 0.002 ± 0.005 


0.24 ±0.01 ±0.02 


0.99 ± 0.03 ± 0.09 


1.04 ± 0.03 




3.3 


- 4.1 


0.031 ±0.001 ±0.003 


0.145 ±0.009 ±0.014 


0.57 ±0.02 ±0.05 


1.06 ±0.03 




4.1 


- 5.3 


0.0157 ±0.0008 ±0.0014 


0.065 ± 0.005 ± 0.007 


0.268 ±0.010 ±0.024 


1.04 ± 0.03 




5.3 


- 7.6 


0.0049 ± 0.0003 ± 0.0005 


0.021 ± 0.002 ± 0.002 


0.095 ± 0.004 ± 0.009 


1.09 ± 0.04 




7.6 


- 12.0 


0.0010 ±0.0001 ±0.0001 


0.0063 ± 0.0010 ± 0.0009 


0.0183 ± 0.0008 ± 0.0020 


1.06 ±0.06 



